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ABSTRACT 

We present images from a long term program (MOJAVE: Monitoring of Jets in AGN with VLBA Experi- 
ments) to survey the structure and evolution of parsec-scale jet phenomena associated with bright radio-loud 
active galaxies in the northern sky. The observations consist of 2424 15 GHz VLBA images of a complete flux- 
density limited sample of 135 AGN above declination -20 degrees, spanning the period 1994 August to 2007 
September These data were acquired as part of the MOJAVE and 2 cm Survey programs, and from the VLBA 
archive. The sample selection criteria are based on multi-epoch parsec-scale (VLBA) flux density, and heavily 
favor highly variable and compact blazars. The sample includes nearly all the most prominent blazars in the 
northern sky, and is well-suited for statistical analysis and comparison with studies at other wavelengths. Our 
multi-epoch and stacked-epoch images show 94% of the sample to have apparent one-sided jet morphologies, 
most likely due to the effects of relativistic beaming. Of the remaining sources, five have two-sided parsec- 
scale jets, and three are effectively unresolved by the VLBA at 15 GHz, with essentially all of the flux density 
contained within a few tenths of a milliarcsecond. 

Subject headings: galaxies : active — galaxies : jets — radio continuum : galaxies — quasars : general — BL 
Lacertae objects ; general — surveys 



1. INTRODUCTION 

High angular resolution studies of the kinematics of AGN 
jets over many years have led to a better understanding of 
the process of acceleration and collimation of the relativis- 
tic jets associated with the ejection of relativistic plasma from 
a central supermassive black hole. Early observations were 



made with ad hoc arrays comprised of relatively few anten- 
nas and focused on a few relatively strong and fast super- 
luminal sources (e.g ., ICohen etaUl 19751: IShaffer etalJll975l: 
iSchilizzi et al.lll9'75l) . In these studies it was difficult to trace 
details of the jet outflow due to limited temporal and inter- 
ferometric coverage of the arrays, which consisted of existing 
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antennas with differing instrumental characteristics. More- 
over, by concentr ating on only the fastest or most inte rest- 
mg sources, fe.g.. lCohen et al.lll977l:IZensus et al.|[T995l) . the 
kinematic results did not reflect the full range of speeds seen 
in the overall population. With the start of VLBA operations 
in the mid 1990s, we were able to obtain a reasonably uniform 
set of multi-epoch 15 GHz (2 cm) observatio ns on 132 AGN 
jKellermann et alj|1998l; IZensus et all 120021) which covered 
the period 1994 to 2002. From these data w e were able to dis- 
cuss the statistical properties of AGN jets ("Kellermann et alj 
|2004; Kovalev et al. 2005; Homan et al. 2006), includin g the 
intrin sic properties of the parent population (ICohen et al.l 
l2007h and the detailed nature of a few source s of particular in- 
terest (iHoman et al. 2003; Vermeulen et al. 2 0031: lOster et alj 
l200irKadler et al..i2008uKovalev et al.>i200"£ 

Although these studies included most known "flat spec- 
trum" (a > -0.5, where Si, cx i/") radio sources with to- 
tal flux density above well defined limits, they included or 
excluded some sources because of the contribution of ex- 
tended emission, errors in the measured spectral index and/or 
spectral curvature. As described in paper I of this series 
dLister & HomanI I2005h . our new MOJAVE (Monitoring of 
Jets in Active Galactic Nuclei with VLBA Experiments) sam- 
ple includes all sources with measured total VLBA flux den- 
sity at 15 GHz greater than a defined flux density at any 
epoch during the period 1994.0 to 2003.0. For the 96 so urces 
in comm on with the studies of iKellermann et all (11998 ) and 
IZensus et al. (2002), the time scale for determining the speed 
and possible accelerations or decelerations is now extended 
from 7 to 13 years. Lister & Hpman (2005) give more de- 
tails of the sample definition and also discuss the single-epoch 
linear polarization prop erties of the MOJAVE sample, while 
iHoman & Listed (l2006h discuss the circular polarization prop- 
erties (Paper II). 

The MOJAVE/2 cm Survey program is among the largest 
multi-epoch AGN VLBI surveys carried out to date, and 
is complemented by several programs at other wavelengths. 
Thes e include the 5 GHz Calte ch-Jodrell Flat-Spectrum Sur- 
vey ( IBritzen et al.ll2007l j^08), the 2 & 8 GHz Radio Ref- 
erence Frame program dPiner et al.1 l2007h . and th e 43 GHz 
Boston University AGN monitoring program ( Jorstad et al.l 
12007) . Each of these programs offers various trade-offs in 
terms of sample size, angular resolution, and image sensitiv- 
ity. Shorter wavelength studies, such as the Boston University 
program, provide better angular resolution than at 15 GHz, 
but the sensitivity is poorer, jet features fade faster, and the 
impact of tropospheric conditions is more severe. Thus, short 
wavelength observations need to be more frequent, and con- 
sequently must focus on much smaller samples of extremely 
bright sources. Longer wavelength studies have better surface 
brightness sensitivity and so can trace the jet motions out to 
greater distances, but at the expense of degraded angular res- 
olution. 

This paper is the fifth in a series describing the results of the 
MOJAVE program, with the most recent ones disc ussing the 
kiloparsec-scale jet properties (Cooper et al."2007')(Paper III) 
and parent luminosity function (Cara& Lister 2008)(Paper 
IV) of the sample. Here we present a complete set of 15 GHz 
VLBA images of the flux-density limited MOJAVE sample of 
135 compact extragalactic radio jets obtained during the pe- 
riod 1994 August to 2007 September The overall layout of 
the paper is as follows: in §|2]we discuss the sample defini- 
tion and selection criteria, in §[3] we discuss our observational 
program and data reduction method, and in § |4] we summa- 



rize the overall parsec-scale morphological properties of the 
AGN jets and describe our current monitoring program. In 
subsequent papers we will discuss their kinematic properties, 
overall demographics, and polarization properties. 

2. SAMPLE DEFINITION 

In Paper I we described our selection criteria for the MO- 
JAVE sample, which is based on compact radio flux density in 
order to provide the best comparison with Monte Carlo sam- 
ples o f relativistically beamed jets (e.g., iLister & Marscherl 
Il997h . By using the milliarcsecond scale (VLBA) 15 GHz 
flux density rather than total flux density as the selection cri- 
terion, we were effectively able to exclude any contribution 
from large-scale emission, leaving a sample consisting al- 
most entirely of radio loud AGN with relativistic jets pointed 
close to the line of sight (i.e., blazars). The only exceptions 
were a handful of nearby radio galaxies and peaked-spectrum 
sources whose jet axes probably lie much closer to the plane 
of the sky. 

The MOJAVE selection c riteria resemble those of the orig- 
inal VLBA 2 cm Survey of IKellermann et all (Il998h . and are 
defined as follows: 

• J2000 declination > -20° 

• Galactic latitude \b\ > 2.5° 

• Total 15 GHz VLBA flux density of at least 1 .5 Jy (> 2 
Jy for sources below the celestial equator) at any epoch 
during the period 1994.0-2004.0. 

Because of the highly variable nature of strong, compact 
AGN, we did not limit the flux density criterion to a single 
fixed epoch. Doing so would have excluded many highly vari- 
able sources from the sample, and would subsequently reduce 
the robustness of statistical tests on source properties. Includ- 
ing the variable AGN also provides a more complete sample 
for comparisons wi th AGN surveys a t other wavelengths, such 
as the 3rd EGRET dHartman et al.ll999.) and upcoming Fermi 
dThompsonll2004l) gaimna-ray catalogs. The higher flux den- 
sity limit for the southern sources was chosen since the VLBA 
has reduced hour angle coverage for this region of the sky and 
consequently poorer imaging capability. 

We note that for several sources in our sample, we did not 
have any 15 GHz VLBA measurements that formally met 
our flux density criterion. However, we were able to infer 
their VLBA flux density at other epochs by using an exten- 
sive database of flux density measurements from the UM- 
RAO and R ATAN tele scopes spann ing 1994.0 to 2004.0 (e.g., 
iKovalev et a l. 1999; Aller. AUer. & H ughes 1992). An impor- 
tant step in this method was to determine the amount of ex- 
tended (non-VLBA) flux density in each source using near- 
simultaneous UMRAOA'LBA measurements made at multi- 
ple epochs. 

In Paper I we originally identified a total of 133 AGN 
that satisfied these criteria. In processing additional archival 
VLBA epochs for the current paper, we have subsequently 
identified two more sources (0838H-133 and 1807H-698) which 
met our selection criteria, bringing the total number up to 
135 AGN. One additional source, 1150H-497, appears to have 
possibly met the MOJAVE criteria based on interpolation of 
single-dish RATAN radio observations in 2003 August, how- 
ever, we did not obtain any 15 GHz VLBA epochs on this 
source prior to 2008. This source is no w currently in our ex- 
tended monitoring program (see § l4.3l l. We believe the over- 
all completeness of the MOJAVE survey to be high, with few 
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sources missed, because of our large pre-selection candidate 
list (see Paper I) and the extensive flux density database that 
was available from the VLBA, RATAN, and UMRAO facili- 
ties. 

We list the general properties of the MOJAVE sources in 
Table 1. The redshifts are gathered from the literature, and 
are currently only 93% complete, mainly because of the fea- 
tureless spectra of several BL Lacs in t he sample. We list the 
nature of the optical counterpart from I Veron-Cettv & VeronI 
([200^ for each AGN in column 6, except as noted in Ta- 
ble 1. At the present time, only four sources (0446H-113, 
0648-165, 1213-172, and 2021H-317) lack pubUshed optical 
counterparts. 

The MOJAVE sample is dominated by quasars, with 
the weak-Uned BL Lacertae objects and radio-galaxies 
making up 16% an d 6% of the sample, respectively. 
IVeron-Cetty & VeronI (|2006') have recently reclassified as 
quasars some BL Lac objects that have occasionally displayed 
emission lin es slightly wider tha n the nominal 5 A equivalent 
width limit (IStickel et al.lll99? ). Since these are still much 
narrower than those typically found in quasars, we have cho- 
sen to list the previous BL Lac classifications for these ob- 
jects, and have indicated these with flags in Table 1. The clas- 
sifications o f the radio spectra l sha pe given in column 7 are 
described bv' Kellermanni et all (12004) . The MOJAVE sample 
is heavily dominated by blazars with flat overall radio spectra, 
which we define as a spectral index a flatter than -0.5 at any 
frequency between 0.6 and 22 GHz. The five steep-spectrum 
sources in the sample have strong extended emission on arc- 
second scales that dominates the integrated spectrum, but their 
parsec-scale emission still met our selection criteria. Two of 
the latter AGN (1458H-718 and 1823H-568) are frequently re- 
ferred to as compact steep spectrum sources in the literature. 
All of the sources originally described in Paper I as peaked- 
spectrum, with the exception of 0742H-103, have subsequently 
turned out to have variable radio spectra and are therefore 
no longer considered a s stable peaked-spectrum sources (e.g., 
iTorniainen etani2005l) . 

3. OBSERVATIONS AND DATA REDUCTION 

The 15 GHz VLBA data presented in this paper consist 
of observations carried out between 1994 and 2007 as part 
of the VLBA 2 cm Survey and MOJAVE programs, and 
data from the online NRAO data archive' spanning the pe- 
riod 1998 - 2006. The multi-epoch observations for each 
source, along with the corresponding image parameters, are 
de scribed in Table 2. The 2 cm Survey observations analyzed 
by iKellermann et al.l (1200 4) and E. Ros et al. (in prepara- 
tion) consist of approximately one hour integrations on each 
source, broken up into approximately five minute-long scans 
separated in hour angle to improve {u,v) plane coverage. A 
similar observing method and integration times were used in 
the MOJAVE observations from 2002 May to 2007 September 
(VLBA codes BLl 1 1, BL123, BL137, and BL149; as flagged 
in Table 2), and are described in Paper I. During 2006 (VLBA 
codes BLl 37), the 15 GHz integration times were shortened 
by a factor of ~ 3 to accommodate interleaved scans at three 
other observing frequencies (8.1, 8.4, 12.1 GHz). The MO- 
JAVE and 2 cm Survey observations were recorded at a data 
rate of 128 Mbps, which was increased to 256 Mbps in the 
epoch 2007 July 3 and epochs thereafter. Beginning with the 



2007 January 6 epoch, the number of sources observed in each 
24 hour MOJAVE session was increased from 1 8 to 25 sources 
to accommodate our larger monitoring sample (see § l4.3l l. 

In order to obtain the best possible temporal coverage to 
study jet kinematics, we re-analyzed all past VLBA 15 GHz 
observations that were available in the VLBA archive from 
1998-2006 with good interferometric data on the sources in 
our sample. Since we require good {u,\') coverage to track jet 
motions accurately, we only considered observations in which 
there were at least three scans on the source, suitably sepa- 
rated in hour angle. We also omitted any archival observa- 
tions that had fewer than 8 VLBA antennas. Any observa- 
tions which lacked data from either of the outermost antennas 
(Mauna Kea and St. Croix) were also omitted if they had 
fewer than 6 scans on the remaining antennas. Finally, we 
excluded a few archival epochs on some sources after pro- 
cessing them if their overall data quality was very poor due to 
weather and/or technical problems. Further details on these 
archival epochs can be obtained from the NRAO data archive 
using the appropriate VLBA code and observation date listed 
in Table 2. 

Including the MOJAVE and 2 cm Survey epochs, we ana- 
lyzed a total of 416 VLBA experiments, yielding 2424 indi- 
vidual source epochs on the MOJAVE sample. A portion of 
these were re-analyzed and/or presented in earlier 2 cm Sur- 
vey and MOJAVE papers. We present the entire set of source 
epochs here for completeness. 

The temporal coverage per source varies between 5 and 89 
epochs, with a median of 15 epochs (column 5 of Table 1). 
The wide dispersion reflects the fact that some AGN are ob- 
served much more frequently with the VLBA, either as cali- 
brators, or because they are sources of particular astrophys- 
ical interest. Also, attempts are made in the 2 cm Survey 
and MOJAVE programs to observe the sources with the most 
rapid motions more frequently, given the constraints associ- 
ated with VLBA's dynamical scheduling system. In setting 
the observing intervals, we have aimed to not let the jet fea- 
tures move more than about half of their typical angular spac- 
ing between epochs. This is done to facilitate the unambigu- 
ous identification of the same features from one epoch to the 
next. 



3.1. Data Reduction 

We processed the data according to the standard proce- 
dures described in the AIPS cookbook". Our overall reduction 
method was to initially flag out bad data on the basis of ex- 
cessively high system temperature (typically due to inclement 
weather), remove sampler bias with the task ACCOR, cali- 
brate the correlator output using antenna system temperature 
and gain curves, and apply an atmospheric opacity correc- 
tion with APCAL. The pulse calibration signals were used 
in all the experiments to align the phases across the IFs. In a 
few cases where no pulse calibration signal was available, we 
performed a fringe fit of a single scan on a bright calibrator 
and used the result to align the delays for the other sources. 
In the case of the archival epochs, and for MOJAVE epochs 
observed after 2007 Aug 9, we deemed the standard global 
fringe-fitting step to be unnecessary because of the strong and 
compact nature of all the sources, their well-determined sky 
positions (e.g., from the ICRF, RD V, and VCS surveys, e.g., 
iFev et al.ll2004t iPetrov et al.]|2008h . and the well-determined 
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^ http://www.aips.nrao.edu 
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positions of the VLBA antennas. After a final bandpass cor- 
rection and a short (30 second) phase self calibration us- 
ing a point source model, we exported the data to Difmap 
dShepherd 1997). The imaging and self-calibration process 
in Difmap was largely automated, with human input required 
mainly for data flagging. The procedure was facilitated by 
the use of a mySQL database containing the locations of clean 
boxes for each source, which we determined during the pro- 
cessing of previous epochs. Since the overall extent of the in- 
dividual sources did not significantly vary over the time period 
covered by our survey, the use of pre-determined clean boxes 
provided additional a-priori information that guided the self- 
calibration process and facilitated rapid convergence toward 
the final model. For experiments in which polarimetric cross- 
hand data were recorded, we performed additional calibration 
steps to correct for antenna leakage terms and the absolute 
electric vector position angle on the sky, as described in paper 
I. These multi-epoch polarization data will be presented in a 
separate paper 

3.2. Images 

We produced the final images using Difmap with a 
multi-resolution clean algorithm developed by Moellenbroc^ 
((1998). The technique involves 'super-resolving' the bright- 
est and most compact features in the initial cleaning stages, 
which leads to a more accurate representation of lower sur- 
face brightness features in the final (natural-weighted) clean- 
ing stages. The final step consisted of a large iteration deep 
clean, which also helps bring out these features. This pro- 
cedure does, however, alter the noise statistics in the final 
cleaned image. For this reason we obtained the rms image 
noise values from a blank sky region a full arcsecond away 
from the CLEAN image center after the final cleaning stage. 
This provided a better estimate of the true rms image noise as- 
sociated with the interferometric visibilities, since this region 
was located well outside the deep CLEAN box of the main 
image. 

We present the un-tapered, naturally-weighted images of all 
the individual source epochs in Figure Set 1 . The dimensions 
of the restoring beam in each image vary depending on the 
(u,v) coverage of the individual observation. Information on 
the restoring beam can be found in Table 2, along with the 
lowest contour level, VLBA experiment code, and observa- 
tion date. The image dynamic range varies with (u,v) cover- 
age and integration time, with typical values of approximately 
5000: 1 . We obtained the total flux densities by summing all 
emission from pixels in the image above the lowest contour 
level within the region shown in Figure 1. Gaussian model 
fits to the visibility data, as well as a kinematic analysis of the 
sample, will be presented in a separate paper 

In order to show the fullest possible extent of the jets, 
we have also produced 'stacked' naturally- and uniformly- 
weighted images of each source (Figure Set 2). We created 
the former by first restoring the image at each epoch with an 
identical beam whose dimensions corresponded to the median 
naturally-weighted beam of all epochs on that source. We 
then shifted the convolved image at each epoch to align the 
fitted positions of the core features, and then combined them 
all with equal weight to produce an averaged (stacked) im- 
age. We created the uniformly-weighted images in the same 
manner, using a restoring beam that was 80% of the naturally 
weighted beam size in both dimensions. 

We note that in the case of a few sources located within 
~ 1 .5 degrees of the celestial equator, the individual baseline 



tracks in the (u,v} plane are almost exclusively in the east-west 
direction, resulting in relatively poorer interferometric cover- 
age. This can create artifacts to the direct north and south 
of the brightest feature, as seen in some of the stacked and 
individual epoch images (e.g., 1055H-018, 2134H-004). VLBA 
images of equatorial sources with jets aligned in the north and 
south directions on the sky suffer particularly from this effect. 

A very small number of our images may also suffer from 
the effects of rapid intraday variability (IDV). Since the imag- 
ing procedure assumes the source flux density is not signif- 
icantly changing on the (several hour long) timescale of the 
observation, intrinsic source variability can add error to the 
gain solutions, with a resulting increase in image noise. We 
ha ve discovered one such case i n 1 156H-295, epoch 2007 Feb 
5 dSavolainen & Kovalevll2008h . We are currently analyzing 
our whole dataset for other instances of IDV (Kuchibhotla et 
al., in preparation). 

The full set of all images are included in the online version 
of this article. Linearly interpolated movies which show the 
jet motions and structural variability, FITS images, and cali- 
brated visibility data are also publicly available on the MO- 
JAVE website . The images presented in this paper cover 
epochs up to the end of 2007. More recent images and cal- 
ibrated visibility data from our current program are uploaded 
to the web site as soon as they are reduced, typically 3 to 4 
weeks after each observing session. All of the raw visibil- 
ity data associated with the MOJAVE program are publicly 
available immediately after correlation from the NRAO data 
archive. 

4. MORPHOLOGY 

The appearance of blazar jets depends not only on their 
intrinsic structure, but also on the Doppler factor associated 
with the highly relativistic motion of individual features in the 
jet. Their small viewing angles can highly exaggerate mild 
curvatures in the jet flow, and relativistic beaming effects can 
amplify or diminish the apparent brightness of different jet re- 
gions as the result of even minor changes in the direction of 
motion with respect to the line of sight. It is thought that all 
jets must be intrinsically double-sided to explain the transport 
of energy from the central engine to the double radio lobes 
commonly found in FRI and FRII radio galaxies and in some 
quasars. In many cases, however, the visible morphology of 
these jets is highly affected by differential Doppler boosting. 
The apparent flux density ratio of the approaching and reced- 
ing jets may reach as high as 10^ in the fastest known jets, 
giving the illusion of a one-sided jet in VLBI images made 
with limited sensitivity and dynamic range. 

All but eight MOJAVE sources in our sample have the 
asymmetric characteristics of one-sided jets, which is in- 
dicative of the high degree of relativistic beaming in the 
sample. Five sources have two-sided jets, of which four 
are low-luminosity AGN (§ |4.2| l. A small number of such 
weakly-beamed sources are expected to be present in our 
flux-density limited sample, by virtue of their proximity. 
We note that earlier authors have speculated on the two- 
sided nature of the p arsec-scale jet in th e luminous ra- 
dio galaxy 2021H-614 (iTschager et al.1 120001: iPeck & Taylo^ 
,2000) . However, since the southernmost feature is more vari- 
able, more compact dL ister 2001), and has a flat spectral index 
jBarvainis & Lonsdale* 199 8) as compared to the other jet fea- 
tures, we consider it as the core feature in a one-sided jet. 

^ http://www.physics.purdue.edu/MOJAVE 
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Fig. 1. — Naturally-weighted 15 GHz total intensity VLBA contour images of individual epoch observations of the MOJAVE AGN sample. The contours 
are in successive powers of \/2 times the base contour level given in Table 2. Because of self-calibration, in some cases the origin may be coincident with the 
brightest feature in the image, rather than the putative core component Usted in Table 3. Clockwise from the top left: Image of 0003-066, 0333+321, 0838-H133, 
1726+455, 1308+326, 0552+398. Note: this is a figure stub, the fuU figure set is available in the electronic version. 
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0007+106: stack of 15 M0JAVE/2cm Nat.Wgt.Epochs 
Bm:1.17x0.53 mas at -6°; Cntrs: V2 x 0.2 mjy/bm 



Relative R.A. (mas) 



0007+105: stack of 15 M0JAVE/2cm Unif.Wgt. Epochs 
Bm:0.94x0.42 mas at -6°; Cntrs: V2 x 0.6 mJy/bm 
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B 6 4 2 -2 8 6 4 2 -2 
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Fig. 2. — Total intensity 15 GHz VLBA stacked-epoch contour images of the MOJAVE AGN sample. The left hand panels represent the unweighted averages 
of the naturally weighted images presented in Figure Set 1, after shifting each map to make the core component coincident with the origin. The right hand panels 
contain stacked uniformly-weighted images with beam dimensions 80% of the naturally weighted images. From top to bottom are stacked images of 0003-066, 
0007+106, and 0016+731. Note: this is a figure stub, the full figure set is available in the electronic version. 
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4.1. Unresolved and Nearly Unresolved Sources 

We have found that three sources (0235+164, 1741-038, 
and 1751+288) are effectively unresolved by the VLBA at 15 
GHz at all epochs, with essentially all of the flux density con- 
tained within a few tenths of a milliarcsecond. Model fitting 
of the visibiUty data for 0235+164 and 1741-038 indicated 
possible structure on a scale of 0.3 mas; however, we could 
not reliably cross identify any components across multiple 
epochs for any of these sources. Another source, 1324+224, 
is mostly very compact, but has a very weak feature located 
approximately 3 mas to the southwest of the bright core fea- 
ture. However, its measured position was erratic across multi- 
ple epochs and we were unable to define any meaningful mo- 
tions in this source. The stacked-epoch image of 1739+522 
contains a weak diffuse structure to the north of a prominent 
compact core, which may be part of a highly curved jet. 

These sources could appear very compact either because 
their jets are atypically weak, or they possess a highly beamed 
core that swamps the jet emission in our limited sensitiv- 
ity/dynamic range images. It is also possible that they may 
be experiencing a period of inactivity in which new jet com- 
ponents are not being ejected. We have witnessed the latter in 
1308+326, which had very little discernible jet structure prior 
to 1997, and now has a pronounced mas-scale jet (Figure Set 
1). 

With the exception of 0235+164, none of the five com- 
pact sources mentioned above show signific ant extended ra- 
dio structure in deep VLA images (Ulvestad et al.l I198U 
iPerlev etal.lll982t ICooper et al.ll2007h . They are flierefore 
well-suited as calibration sources for flux density and abso- 
lute polarization position angle calibration when simultane- 
ous single-dish/VLA and VLBA data are available. Since any 
observable jet structure in these objects must be on angular 
scales just at or less than the size of the 2 cm VLBA restor- 
ing beam, their kinematics would be best studied with shorter 
wavelength ground and/or space VLBI, which can provide 
better spatial resolution. 

4.2. Two-sided Jets 

Multi-epoch observations of two-sided AGN jets are ex- 
tremely useful for testing whether apparent motions of jet fea- 
tures do indeed reflect the underlying flow. In the case of ran- 
dom pattern speeds, we would not necessarily expect motions 
in the jet and counter-jet to be correlate d. Detailed kinematic 
studies of two-sided AGN jets (e.g., iTavlor & Vermeulenl 
[1997; Vermeulen et al. 2003) have revealed similar speeds in 
the jet and counter-jet, and have provided additional estimates 
of the viewing angles to these jets. In this section we provide 
brief descriptions of the five two-sided AGN jets in the MO- 
JAVE sample. 

0238-084 (NGC 1052): The morphology and kinematics of 
this source have been studied extensively by I Vermeulen et af] 
(12003 ). The jet shows numerous features moving outward 
along the jet and counter-jet at s peeds up to ~ 0.3 c. Multi- 
frequency VLBA ob servations bv' Vermeulen et afl (1200 3) and 
iKadler et al.l (120041) indicate that the core feature is obscured 
at 15 GHz by strong free-free absorption associated with a 
circumnuclear torus. A more detailed analysis of the multi- 
epoch structure of NGC 1052 based on MOJAVE and addi- 
tional 22 and 43 GHz data from a dedicated multi-wavelength 
monitoring campaign is currently underway (E. Ros et al., in 
preparation). 

0316+413 (3C 84, NGC 1275): Although pre- VLBA ob- 



servations indicated a one-sided jet, albeit with complex struc- 
ture, the earliest 15 GHz VLBA observations revealed a previ- 
ously undetected component which had been obscured at the 
lower frequencies u sed previously. Extensive observations by 
IWalker et al.l (12000*) have indicated that the northern lobe is 
likely being affected by free-free absorption from ionized gas 
associated with an accretion disk. Th e central feature con tains 
complex sub-structure and motions dPhawan et al.lll998i) . in- 
cluding several distin ct circularly polarize d regions seen in 
the MOJAVE data (Homan & Listed 120061) and at lower fre- 
quencies (Homan & Wardle 2004). 

1228+126 (M87): The 15 GHz VLBA kinemat ics of this 
nearby AGN were reported by Kovalev et al.l (12007) . and con- 
sist of relatively slow speeds below 0.07c. We have taken 
the core to be coincident with the brightest component, for 
which we have also detected circular polarization (Pap er II); 
a faint counter-jet is observed not only at 2 cm (Kovale v et al. 
, 2007) but also with higher resolution at 13 and 7 mm ( Ly et al. 

1413+135: This highly unusual two-sided BL Lac object 
with an apparent spiral h ost galaxy has been the su bject of 
previous VLBI studies by iPerlman etaP (fT99l [T996l " l200l 
and Wiik et al. (2001). At 15 GHz we see a mostly one-sided 
jet pointing toward the west, but at longer wavelengths there is 
a prominent but more diffuse jet pointing toward the east. The 
line of sig ht to the nucleus shows HI and molecular absorption 
( IConwav|[T999i) . It is the only high-luminosity two-sided jet 
in the MOJAVE sample. Although 1413+135 has been classi- 
fied as a compact symmetric object by Gugliucci et al. (2005), 
its core component is highly variable, which is not a general 
characteristic of that class. 

1957+405 (Cygnus A): The parsec-scale jet of this nearby 
radio galaxy contains a number of well defined features in 
both the main jet and the counter-jet. The jet kinematics have 
been studied by Bach et al. (2005), based in part on 2 cm 
VLBA survey and MOJAVE data, and their multi-frequency 
VLBA analysis suggests that like NGC 1052 and 3C 84, the 
counter-jet is partially obscured by an absorbing torus. 

Although it is tempting to interpret their morphology and 
absorption as evidence that these two-s ided jets are plane- 
of-the sky versions of the one-sided jets, ICohen et alJ (l2007h 
have demonstrated that the low-luminosity two-sided sources 
apparently do not differ from the high lumino sity one-sided 
sources solely because of their orientations. iKovalev et al.l 
(|2007) have also argued that the prominent parsec-scale jet 
in M87 may be intrinsically one-sided. The kinematic inter- 
pretation of the two-sided jets in the MOJAVE sample will be 
presented in a separate paper 

4.3. Current monitoring program and additional sources 

The flux-density limited MOJAVE sample is statistically 
complete on the basis of compact emission at 15 GHz, which 
is essential for understanding the complex selection biases as- 
sociated with observed blazar samples. However, new high 
energy observatories such as Fermi are likely to detect a 
broader range of blazars than those derived from a strictly 
radio-selected sample. For this reason, in 2006 we expanded 
the MOJAVE monitoring program beyond the flux-density 
limited sample to encompass 65 additional sources of inter- 
est. These include i) a subsample of lower-luminosity AGN 
with 8 GHz VLBA flux density greater than 0.35 Jy, 8 GHz 
VLBA luminosity < 10^* W Hz"' and z < 0.3 that we identi- 
fied using the VLBA Calibrator Survey (Bea slev et"ani2002l; 
.Fomalont et al..2003; Petrov et al..2005.) . ii) several additional 
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AGN ident ified as high-confidence EGRET gam ma-ray asso- 
ciations by Sowards-Emmerd et alj (120031 l2004l) . and iii) sev- 
eral AGN known to have particularly interesting or unusual 
kinematics from our pri or 2 cm VLBA survey observations 
dKellermann et a l."2004'). The current AGN monitoring list 
can be found at the MOJAVE website. 

As of Fall 2008, the MOJAVE program is regularly mon- 
itoring 200 AGN with the VLBA at 15 GHz, the UMRAO 
26 m radio telescope at 5, 8, and 15 GHz, and the RATAN 
600m telescope at 1, 2.3, 4.8, 7.7, 11, 22 GHz. We antici- 
pate adding up to 100 additional AGN to our current moni- 
toring list, based on new gamma-ray detections made by the 
Fermi LAT instrument (A. A. Abdo et al., in preparation), 
provided these sources have sufficient compact flux density 
100 mJy) for imaging with the VLBA at 15 GHz. This 
extended sample will ensure complete coverage of several re- 
gions of the gamma-ray/radio flux density plane, and will be 
used for correlative studies. 



5. SUMMARY 

The MOJAVE program represents the largest full- 
polarization high resolution monitoring program of powerful 
AGN jets carried out to date. The complete flux-density lim- 
ited sample of 135 sources contains the brightest and most 
highly beamed AGN visible in the northern sky. Composed 
primarily of blazars, it is especially useful for comparisons 
with studies at other wavelengths, especially at very high en- 
ergies, where blazars tend to dominate the all-sky catalogs. 
In 2006, we expanded the monitoring sample to include an 
additional 65 AGN on the basis of their known gamma-ray 
emission, low luminosity, and/or unusual jet kinematics. This 
list will be further supplemented with new AGN detections 
from the Fermi gamma-ray observatory. 

By combining archival data from the 2 cm Survey and other 
VLBA programs, we have obtained 2424 individual 15 GHz 
VLBA images on the flux-density limited MOJAVE sample, 
covering the period 1994 August to 2007 September The 
full set of images is publicly available on the MOJAVE web- 
site. These data indicate that 94% of the jets are one-sided, 
likely due to relativistic beaming effects. Of the remainder. 
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TABLE 1 

The Complete Flux-Limited MOJAVE AGN Sample 



Source 


Alias 


R.A. 


Declination 


N 


Opt. CI. 


Radio 


z 


Reference for z 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 



0003-066 
0007+106 
0016+731 
0048-097 
0059+581 
0106+013 
0109+224 
0119+115 
0133+476 
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MOJAVE. V. Multi-epoch VLBA Images 



TABLE 1 

The Complete Flux-Limited MOJAVE AGN Sample 
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TABLE 1 

The Complete Flux-Limited MOJAVE AGN Sample 



Note. — Columns are as follows: (1) lAU Name (B1950.0); (2) alternate name; (3) right ascension (J2000); (4) declination (J2000);(5) number of 15 GHz VLBA epochs 
analyzed; (6) optical classification according to theiVeron-Cettv & Veron 1 2006) catalog (with exceptions as noted), where Q = quasar, B = BL Lac object, G = active galaxy, and U 
= unidentified; (7) description of radio spectrum fromt^Kellermann e t alj (2004) , where F=flat, S=steep, and P=peaked; (8) redshift; (9) literature reference for redshift. 

^ Source classified as a quasar in the lVeron-Cettv &~Veron 1 2006) catalog. 

^ Source classified as a possible BL Lac object in the lVeron-Cettv & VeronI 120061) catalog. 

Source classified as a galaxy in the lVeron-Cettv & VeronI 120061) catalog. 

^ Source not listed in the lVeron-Cetty & Verot] 120061) catalog. 

^ Compact steep spectrum source. 



TABLE 2 

Summary of 15 GHz Image Parameters 
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^maj 


-^min 


Bpa 


Itot 


rms 


■l-base 


Fig. 


Source 


Epoch 


Code 


(GHz) 


(mas) 


(mas) 


(°) 


(Jy) 


(mJy bm"' ) 


(mJy bm~^ ) 


Num. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 




0003-066 


1995 Jul 28 


BZ014'' 


15.3 


1.17 


0.51 





2.172 


0.5 


3.0 


1.1 




1996 Oct 27 


BK37D'' 


15.4 


1.20 


0.45 


-7 


1.893 


0.4 


0.6 


1.2 




1998 Oct 30 


BK52D'' 


15.4 


1.36 


0.44 


-14 


2.642 


0.4 


1.5 


1.3 




1999 Jun 17 


BG091 


15.4 


1.36 


0.53 


-9 


2.707 


0.6 


1.4 


1.4 




2000 Jan 1 1 


BK68D'' 


15.4 


1.83 


0.50 


-14 


2.431 


0.3 


4.4 


1.5 




2001 Jan 21 


BK68H'' 


15.4 


1.62 


0.58 


-12 


1.878 


0.3 


0.9 


1.6 




2001 Oct 31 


BR77A'' 


15.4 


1.59 


0.52 


-12 


1.559 


0.2 


1.4 


1.7 




2002 May 19 


BG121A 


15.3 


1.47 


0.57 


-8 


2.192 


0.5 


1.2 


1.8 




2003 Feb 5 


BL111E-' 


15.4 


1.32 


0.53 


-6 


2.847 


0.2 


0.7 


1.9 




2004 Mar 22 


BG144A 


15.4 


1.34 


0.54 


-3 


3.224 


0.5 


1.8 


1.10 




2004 Jun 1 1 


BLlllM" 


15.4 


1.30 


0.49 


-7 


3.302 


0.3 


1.5 


1.11 




2005 Mar 23 


BL123D'' 


15.4 


1.33 


0.53 


-5 


3.034 


0.3 


0.9 


1.12 




2005 Jun 3 


BL123H'' 


15.4 


1.37 


0.53 


-8 


3.017 


0.2 


0.9 


1.13 




2005 Sep 16 


BL123M^ 


15.4 


1.33 


0.52 


-6 


2.718 


0.2 


0.5 


1.14 



Note. — Columns are as follows: (1) lAU Name (B 1950.0); (2) date of VLBA observation; (3) VLBA experiment code; (4) observing frequency in GHz; (5) FWHM major axis 
of restoring beam (milliarcseconds); (6) FWHM minor axis of restoring beam (milliarcseconds); (7) position angle of major axis of restoring beam; (8) total I flux density (Jy); (9) 
rms noise level of image (mJy per beam); (10) lowest I contour (mJy per beam); (11) figure number. 

^ Full polarization MOJAVE VLBA epoch 

^ 2 cm VLBA Survey epoch 



